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Abstract 

In this article, an analytical description is presented based on the master equation. 

This master equation is formed from the system of superconducting phase qubit which 
is coupled to a torsional resonator and damped by a dispersive reservoir. An analytical 
approach for searching of some physical phenomena on the qubit system is presented, 
such as qubit inversion, purity and negativity. In addition, these phenomena are dis¬ 
cussed for the resonance and off-resonance cases for many different initial states. From 
computational results, it is found that the mentioned phenomena depend on the disper¬ 
sive reservoir parameter which leads to their death. However, the complete destruction of 
system coherence in the presence of dispersive reservoir leads to the death of entanglement 
between the qubit and torsional resonator. 

keywords: Superconducting Phase Qubit; Density Matrix; Entanglement; Dispersive Reser¬ 
voir. 


1 Introduction 

In the recent years, the quantum entanglement and mixedness are two biggest challenges 
of quantum information theory. Quantum entanglement is a very important quantum char¬ 
acteristic that has no classical counterpart. Also, quantum entanglement is necessary for the 
realization of quantum communication and the most important computational tasks, such as 
quantum teleportation [^, superdense coding and so on. In addition, it is crucial to under¬ 
stand the fundamental aspects of quantum mechanics, including quantum measurement and 
quantum decoherence [^. Consequently, no system can be completely isolated from its envi¬ 
ronment. Thus in real experimental situations due to the coupling environment, the entangled 
state inevitably loses its purity and becomes mixed. 

Decoherence phenomenon is the most dangerous obstacle for all entanglement manipula¬ 
tions. However, the recent studies of decoherence effect for different systems are increased due 
to the variety of applications in areas as quantum optics and quantum information In 

fact, the interactions between a quantum system with the environment are caused two types 
of irreversible effects, e.g., quantum dissipation and quantum decoherence [^[^. The former 
will cause energy dissipation, and the latter will make the system degenerate from coherent 
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state to classical state. Entanglement evolntion under decoherence has attracted much inter¬ 
est and many researchers have studied it extensively based on various models in the view of 
environment-induced decoherence fTolITT . 


In the previous years, some researchers have found that state entanglement can be decreased 
to zero abruptly and remained zero for a finite time. This phenomenon is called entanglement 
sudden death 12 -16 , where entanglement of two-qubit decays to zero in a finite time under 


the influence of pure vacuum noise [1^ . Precisely, entanglement sudden death has been experi¬ 
mentally observed in an implementation using twin photons 18 and atomic ensembles 19 . In 


many efforts have been devoted to study various types of superconducting (SC) qubit which is 
considered as good candidate to solve this problem. 

Nowadays, superconductors are viable elements and feasible within current experimental 
technology. SC qubit is represented a unique and its interesting approach to quantum infor¬ 
mation and quantum computation because it naturally allow strong coupling. Some aspects 
which are discussed through these systems: entanglement between a flux qubit and a super¬ 
conducting quantum interference device (SQUID) have been realized experimentally. In 
addition, there are other aspects: the strong coupling of a single photon to a SC qubit has been 


experimentally demonstrated 21 by using a one-dimensional transmission line resonator 22 


Recently, some quantum logic gates have proposed 23 by using a system consisting of a num¬ 
ber of SC-charge qubits coupled to a resonator to study different quantum operations. Also, 
a new physical scheme for implementing a discrete quantum Fourier transform is proposed via 
SC qubits coupled to a single-mode SC-cavity 24 . However, the effect of environment cannot 
be neglected in these situations [^, because, the qubits are never isolated, but under some 
environmental influence. 

Our purpose of this article is to give an analysis of the entanglement dynamics of a SC 
phase qubit coupled to a torsional resonator. Consequently, some phenomena, such as the 
qubit inversion, purity and negativity are theoretically investigated, mainly because such states 
are more physical and experimentally relevant. Thus we take into account decoherence and 
analyze its effect on these phenomena. Despite the complexity of the problem, the analytical 
descriptions of the exact solution of the master equation is presented. 

This article is arranged as follows: The physical system of qubit with torsional resonator and 
its dynamics is provided in section 2. The influence of phase-damping on a SC qubit inversion 
and measures of entanglement is presented in section 3. Finally,the conclusions is demonstrated 
in section 4. 


2 The physical system of qubit with torsional resonator 
and its dynamics 

We introduce a theoretical description of the superconducting phase qubit coupled to a 
torsional resonator. A phase qubit with k excess Cooper-pair charges consists of two super¬ 
conducting islands connected to a superconducting electrode through two identical Josephson 
tunnel junctions of a small size with the same capacitance Cj, see Figj^ If the phase qubit is 
working in a regime with <C Eq <C Ej, where Ej, Ec, ks and T are Josephson energy, 
charging energy, Boltzmann constant and temperature, respectively. The Hamiltonian of the 
superconducting phase qubit can be described by the form 

where 0, /, /q, and fc are the phase difference across the junction, bias current, junction critical 
current, and the external flux which is generated by a classical applied magnetic field, respec- 
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Figure 1: Circuit diagram of the phase qubit; a single Josephson junction with Josephson 
coupling energy Ej and Josephson capacitance Cj is driven with a bias current I. The rectangle 
in this diagram refer to dissipative part of the admittance of the bias current I. 


tively. Moreover, k and 0 are quantum mechanical conjugate variables and satisfy [k, 0] = —i. 

In this study, the bias current I is assumed to be zero. The Josephson energy, and the 
charging energy for each junction also are assumed to be Ej = h/o/2e ~ 6.8 x 10“^^ J and 
Ec = 8.3 X 10“^^ J as estimated in (271. Therefore, the effective Josephson coupling 

Ej^ = Ej cos(7r/c) of phase qubit is controlleoBy the external flux fc- If the phase different 0 
is confined to the period tt/c, the states of a single qubit are restricted to a two lowest-energy 
states l^f) and |e). Hence, the Hamiltonian Q can be converted to a spin representation as 
follows: 


Hqu 




-CFz 


( 2 ) 


where Uz 


^ 2 EcEj^ frequency of the phase qubit, and is the usual Pauli spin 

operator. 

The torsional resonator mode is described by a harmonic oscillator-like Hamiltonian as 
follows (28l 


77, 


p2 1 

— -E - T , , 2 t 2 

mode r) T ' ’ 

"-7 777 , ^ 


( 3 ) 


where Py is the angular momentum conjugate to T and T G [0,27r]. The fluctuations of the 
angle T can be described by the fluctuation Tq = y/h/ZnAd in the ground state, where 7^, 
and ojy are the rotational moment of inertia of the torsional resonator, and vibration frequency, 
respectively. 

If we put the SC qubit in the torsional resonator, the effective flux is /c = / sin T. Hence, T 
is measured relative to the direction of /c, Ej is the largest energy scale and coupling between 
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the qubit and torsional vibrational mode can be a large quantity. Thus, we consider the qubit 
is coupled to resonator mode. Also, the external magnetic field is parallely applied to the phase 
qubit plane. Then the external field in this case is given by: /c = / (T — Te), where is 
the angle at equilibrium measured from the direction of fc- In this work, / is the maximum 
value of magnetic fiux when the field is perpendicular to the qubit plane. If Te is non-zero, the 
coupling strength is decreased by factor sinTe, so consider the angle Te = 0. However, one 
concludes that the total Hamiltonian restricted to the two states is given by: 

Ht = - -az + G{a + d^)ax, (4) 

where (ol) o and (fil) d are the (creation) annihilation operators for the vibrational mode and 
phase qubit, respectively, G = is fhe constant coupling energy between the phase 

qubit and the torsional resonator, and co = \JojI-\- 27r^f^Ej/Im is the torsional oscillator 
frequency. Because of the negative sign of the second term, the above Hamiltonian Q in the 
counter-rotating wave approximation can be written as follows: 


H 


hjj 


fujjd)d - -dz + G{ad- + d^d^). 


( 5 ) 


Under the previous hypothesis, the non-unitary time evolution of present interaction system 
with taking into account the dissipative effect is given through the density operator. The density 
operator is described by a master equation of the form: 


hp = i[p,H] + ^ (dzpdz - p 


( 6 ) 


where 7 is a positive parameter that represents a dispersive reservoir (or phase damping of 
the qubit system). It is one of the variables that does not interchange energy with the qubit 
system. Dispersive reservoir of the qubit may occur due to dephasing interactions that might 
arise. For instance, the dispersive reservoir may occur from the elastic collisions in a qubit 
vapor. The exact solution for Eq§ in the case of a high-Q torsional resonator ( 7 /G <C 1 ) is 
obtained by analytic method. The used analytic method is presented by the dressed-states 

This representation is consisting of the complete set of Hamiltonian 


representation 129-31 


eigenfunctions. Here, the states \g) and |e) are corresponding to the higher energy and the 
lower energy, respectively. In the invariant subspace spanned by |e,n -|- 1) and \g,n), the 
eigenvalues and eigenfunctions of the Hamiltonian i§ are given by: 




±1 




H\e,0) = 


hujz 


|e, 0 }. 


( 7 ) 


with the abbreviations = hw{n+^)±Rn, |T±) = /3^\g,n)±/3^\e,n+l), {3^ = 

Rn = G\J (n -|- 1 ) -|- ( 7 )^, and 5 = detuning parameter which represents the 

difference between transition energy of qubit and torsional energy. 

To obtain the analytical solution for Eqj^ we suppose that the qubit is initially in the 
general superposition of the excited state and ground state while the torsional resonator is 
assumed to be initially in a coherent state, that is 

p(0) = ^cos6*|e)-|-e*‘^sind|(yf)^ ^cos6*|e)-|-e*‘^sind|(yf)^ 0 |q;)(q;|. (8) 
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Consequently, the time-dependent analytical solution for the density matrix with the initial 
previous state is given by: 


OO 

p{t) = Pee|e,0)(e,0| + +1| 

n,m=0 

+Pe 5 ^^’'"(^)|e,rr + +Pg+^’"*+^(t)|e,n+ l)(e,m + 1||, (9) 


where all elements of the density matrix ([^ are also functions in 6, cp, 5. The elements of the 
previous density matrix are given by the following formulas: 

x„a:^e-*^"'-*{2/3+/3-/3+/5-^sin2 0[E+- cos(i?+-1) - K++ cos(i?++1)] 

Hn,m COS^ + /5T^/3“^e*<-*)E+- 

+(/^n ] 

- /3+^e-*<-*)K++] 

+ Vm+i" - /5+^e*<-*)K++]}, 

V m; 

|x„P{sin^ 0[n+“ - n+“] -h cos^ 9[uij+~ + n+“] -F r2+“ -F r+“}, 

V n=m. 




‘w = 1 




{sin^ d[(/5- 

+(/5++ /5-V+'e"'^"’^')K++] 

+2/3+/3T/5;J;/5“en,m cos2 0[E+- cos(i?+-1) - K++ cos(i?++1)] 


pr (^) = 


<^/^n sin 26 [— -_|_ ] 


I <^C'n C'n r Q—^p 

sin 2^ |-^^_2^ 


2\/?Tl-pT 


_2 — f 


V UT^ m; 

a:„p{sin^ 6*[tr7+’“ -|- n+’“] -|- ^n,n cos^ P[n+’“ — n+’“] — f 2 +’“ — r+’“}, 

V n=m. 


n+l,m 

Reg 


{t) = 


XnXl^e s\v?e[{^^ e 


+ 2 ^ii?+^t)E+^ + (/ 5 +^e - / 3 ^^e*<-*)K++] 


— /5~^ 6' 

n c"n " l\p"m /-^n,m ' VC"m 

+i^[(/5;["/5“'e-'<-' + /5T'/3+"e<-*)E+- + (/5+"/3+'e-<-* 

“ K++cos(i?++t)]}, 

V n^^ m; 

knlM^n sm^ ^[(/^n ^ E+- + - /5T^e*<"*)K++] 

+/^n/^n ^n,n COS^ d[{^R - ^T^)e^E+- + 

I /3+^/3n ^sin26>(ct+Q;*) — 

sin2^ \a ^ P+^^-iRiXt ^ Q-^^iRtt.t 
\/n+l L 2 

V n=m 


e* 

- /5+"e-*<-*)K++] 


')KJJ] 


+^iy/ 3 ^*e-‘''"-“‘ + A;V“”."‘) - a*a^ft 7 ^cos(flJ-+()]K.t-+}, 
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Where 


^^m+l ^ (^n+l,m)t^ Pee = e '"I" COS^ P, a = ae Cn,m = 




Xn = e(-^)^, R+-=R^-R^, R++=-fin + Rm, 

= |[1 - “ /^n“ /^n 

Km = |[i + (/^n + /^n V,;' - K'K' - 

E+- K++ 6 „ = 4/3f/5-^ =/^^(n - m), 

p+- ^ pt^nsin2B — /5+^(Q;e“*^"’”* + Q;*e*^"’”*)]]K++, 

^ _ /5-^)e^E+-, 

n+’" = [{K^ + /5“^) smh(3^) + 2l3fpf cosh(3^)]E+-, 

= [{K + K) cosh(3^) + 2pfl3f sinh(3^)]E+- and 
= 2 ^f f 3 f cos{R+%t)K+%. 


■y/(n+l)(m+l)’ 


To calculate the asymptotic behavior of the density matrix, we set Gt ^ oo, as follows: 


-| LAkJ 

p{Gt^oo) = p“|e,0)(e,0| + -J]^|a:„p(^sin^P + ^„,„cos^P^ 

n=0 

x(|e,n + l)(e,n + 1| + |^,n)(^,n|), 
which is a statistically mixed state. 


3 The influence of dispersive reservoir on some physical 
phenomena 

In the current section, we study the existence of some phenomena of the previous system. In 
an effort to present a numerical characterization, we have performed some calculations for the 
properties quantity for a particular set of parameters, some of which can be considered as 
realistic, while some other parameters look perhaps too optimistic. 

3.1 SC qubit inversion phenomenon 

The qubit inversion is one of the important dynamic variables of the properties of quantum 
mechanics. It gives information about the behavior of the qubit during the interaction period. 
Therefore, by using the analytic solution of the master equation and keeping in mind that the 
state \g) is the one with higher energy, we discuss the effect of dispersive reservoir on qubit 
inversion, which is given by: 


W„(t) = pS + E('’«W-4w)^ ( 10 ) 

i 

It is generally accepted that the revivals of inversion appears as a consequence of quantum 
coherence which are built up during the interaction between the resonator and the qubit. 

Numerical calculations which indicate the influence of dispersive reservoir ^ on qubit in¬ 
version is exhibited in Figj^, Figj^ and Figj^. The time evolution of Wgu is plotted as a 
function of the scaled time — and T, in the resonance case, for a different values of 6 and a. 
In the absence of detuning effect and see FigJ^ and Figl^, one finds that Wgu fluctuates 
around zero, and the collapse and revival phenomena start to appear. It is seen from Figj^ 
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and Figj^ a short period of collapse after onset of the interaction followed immediately by a 
reasonable revival period. After this period of collapse, one can not see another collapse period 
because of existence of interference in the fluctuation, so in these cases of 9 the revival period 
is not apparent. The reason for this is due to the value of coherence parameter a. If a = 6 , the 
revival and collapse phenomena are more pronounced than the other case for different values 
of 0, say |. This result is observed in Figj^. Furthermore, for many different values of 6, the 
revival and collapse phenomena are more pronounced. Also, it is observed from Figj^that the 
periodic oscillations of Wgu between 1 as a maximum value and —1 as a minimum value around 
0 . 

If the dispersive reservoir parameter is regarded § 7 ^ 0, the time evolution of Wgu is shown 
in Figsj^, Figj^ and Figj^. At a weak dispersive reservoir, the collapses and revivals phe¬ 
nomenon is still observed for a short time. The increase of collapses region slowly due to the 
influence of the damping. But with large values of dispersive reservoir, the amplitude of revivals 
is greatly decreased. With ^ > 0.1, all revivals are completely washed out, and the system is 
collapsed completely. This means that the probability of existence for the qubit in the state 
\g) and |e) are equal. 

To investigate the influence of detuning parameter ^ and ^ parameter on the behavior 
of qubit inversion, we display in FigJ^ the behavior of inversion in off-resonance case with 
6 * = |, ^ = 6 and a = 3. Firstly, for ^ = 0, one can observe from Figj^, the revivals and 
collapses phenomena are more pronounced than the case of ^ = 0 , and the curve of qubit 
inversion is shifted down. Also, one can note that from Figj^ the revival period is increased, 
and it can distinguish other periods of revival in that case, that is in the case of off-resonance 
for this case of 9 the revival period is more apparent. In fact the revival time is given by: 
tr = 2ttG which increases by increasing S. However, we can see from Fig, 2 d, that 

the detuning parameter leads to the state of qubit to tend to \g) more than |e). Besides that, 
if ^ 7 ^ 0, the collapses and revivals phenomenon is still observed in a very short time. This 
phenomenon disappears slowly by increasing time due to the influence of One can observe 
rapid deterioration for the revivals of Wgu at the large values of From Figj^, we observe 
increasing in the behavior of Wgu curve, after that the value of Wgu is almost to zero at time 
Gt ~ SOtt. Moreover, the collapse region has becomes very large. However, if we take large 
values of ^ for the different values of 9, the amplitudes of the fluctuations are diminished, the 
collapse periods are elongated, and a large decrease is happened in the amplitudes of periodic 
oscillations of maximum value and minimum value for the inversion. 


3.2 Dynamical properties of coherence and entanglement phenom¬ 
ena 

In order to understand the influence of dispersive reservoir of qubit on some dynamical 
properties, such as coherence and entanglement, we use in the following the solution in Eqj^ 
Then we investigate the effect of dispersive reservoir on the dynamical properties. 

3.2.1 Dynamical properties of coherence phenomenon 

In this part, we study the behavior of purity as a measure of the amount of coherence loss 
caused by the dispersive reservoir of qubit and the unitary qubit-vibrational mode interaction. 
The purity of this state is measured by the von Neumann entropy as follows: 

00 

Squ{r){'t) — ~ ^ ^ ^qu(r) (^) ^qu(r) (^) ^ 

i=l 
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Figure 2 : The evolution of SC qubit inversion plotted against the scaled time ^ and dispersive 
reservoir ^ in the resonance case for 6 * = ^ in (a), 6 * = | in (b) and 6 * = | in (c), and off- 
resonance case for 6* = I and ^ = 6 in (d), with the coherence parameter a = 3 in (a, b, d) 
and a = 6 in (c). 










where are the instantaneous eigenvalues of the two subsystems (qubit and torsional res¬ 

onator). In a bipartite quantum system, the system and subsystem entropies, at any time t, 
are limited by the following Araki-Lieb inequality — Ar] < S < |*S'g„ + 5'r|. As a result, if the 
density operator p describes a pure state, the entropy S becomes 0 and remains constant. This 
means that, if the system is initially prepared in a pure state at any time t > 0, the reduced 
entropies of the two subsystems are identical, that is, Squ{t) = Sr{t). Otherwise A > 0, that is, 
Squ{t) 7 ^ Sr{t). 

Here, we begin from pure state, therefore, we either use Sr or Sgu in fhe absence of dispersive 
reservoir, to measure the amount of entanglement between the two subsystems, whether in the 
resonance or off-resonance cases. Exactly, S',j„(0) = S'j.(O) = 0. In this work, we devote our 
study to the von Neumann entropy for the reduced density operator of the qubit that can be 
expressed as follows: 


Sj.w = (12) 

i=± 

where (f) = + Er=o[Pef (^) + Pg'^i^)]) ± with A„(t) = \(pf^ + 

E^o[dee’^(^) ^n{t) = \ our discussiou for the influence 

of dispersive reservoir of qubit on the purity, the time evolution of the purity against ^ and 

in the resonance case with 6 * = |, (p = 0 and a = 3, is displayed in Figj^. From Figj^, 
if ^ = 0, one can observe that the behavior of curve of Squ is not periodical. In addition, 
we see that the time evolution curve of qubit entropy has local extreme at the revival time 

{tn = = 3 -y/s/mUi /EjiOz f ^ ^, indicating that the qubit state does not return most closely 

to its initial state at this time. Also, one can observe that the value of minima for Squ in Figj^ 
is increased when the time increases. Finally, the Squ oscillate around its maxima as time 
becomes very large. Besides that, in the presence of the dispersive reservoir ^ 7 ^ 0, we cannot 
use the Squ as a measure for entanglement but we use it as a measure to the amount of the 
coherence loss of the qubit state. At a weak dispersive reservoir, one can find the amplitudes of 
the oscillations of entropy diminish rapidly, and the local extreme of the purity at half-revival 
time for the mentioned measure are still apparent, where at this time the qubit and vibration 
mode are approximately disentangled. If ^ = 0.1, the amplitudes of the local extreme of 
the purity oscillations are more rapidly deteriorated. Due to damping, the amplitude of the 
purity of qubit is suppressed. Therefore, the qubit coherence loss may occur due to dephasing 
interactions of qubit with the vibration mode. However, the purity loses its coherence more 
than its gain due to dispersive reservoir, and hnally completely loses coherence and fall into a 
statistically mixed state. 

As well as, we investigate in Figj^ and Fig|^,|^), the evolution of purity in off-resonance 
case with ip = 0 and a = 3 for a different values of 6 and In the absence of see 
Figj^, Figj^ and Figj^, one can note that the detuning parameter affects the fluctuations 
in the temporal evolution for Squ, and results in deterioration of the oscillations. One can 
also note from these figures, the local extreme for the mentioned measure are pronounced and 
still apparent in all curves. However, we see that the amplitudes of oscillations of purity are 
decreasing due to the increase of the value for ^ in these different figures, that is the detuning 
parameter improves the influence of dispersive reservoir on this phenomenon. In addition, when 
^ 7 ^ 0 , one can And from Figj^, Figj^ and Figj^ at a weak value of that the focal extreme 
for the mentioned measure is about to disappear. Also, one can observe at the high values of 
that the amplitudes of oscillations for all curves are deteriorated completely, and the purity is 
increased monotonically and no longer equal zero due to the influence of the dispersive reservoir. 
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Figure 3: The evolution of purity plotted against the scaled time ^ and dispersive reservoir 
in the resonance case for 0 = | in (a), and off-resonance case for d = | and ^ = 4 in (b), 
= I and ^ = 10 in (c) and 0 = ^ and ^ = 25 in (d), with a = 3. 
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however, from Figj^ for 9 = \ (phenomenon of coherent trapping) and ^ = 10 , and from 
Fig,^ for 6 * = I and ^ = 25. We note that the behavior for two curves of the purity in 
these cases of detuning is periodical, and note that the oscillations periods for two curves are 
apparent after a large period of time on the case in which ^ = 4. Then, we conclude that 
the behavior of curves of purity is also periodical for the values more than ^ = 25 for the 
different values of 6. Finally, when ^ = 1, we get A^(G't cxd) = ^ ± cos^ d, then 

SquiGt oo) ~ ln2 = 0.693, which means that the asymptotic qubit state is a statistically 
mixed state, that is the mixedness of the phase qubit state is complete. 


3.2.2 Dynamical properties of entanglement phenomenon 


Because of the existence of dispersive reservoir parameter in Eqj^ the purity cannot be a 
measure for the entanglement. Therefore, we use the negativity N(p) to measure the amount 


of the entanglement that can be expressed as 33 


N(p) = 


max 



(13) 


where A* is a negative eigenvalue of the partially transposed density matrix of the vibrational 
mode-phase qubit. Here, A^ is computed by using numerical calculations. For an entangled 
mixed state, N{p) is positive. If N{p) = 0, the states are separable. Moreover, the negativity 
is used to estimate the amount of entanglement of the final state if we start with a pure or 
mixed state. 

The infiuence of dispersive reservoir of qubit on the negativity as a function of — and a 
parameter | in the resonance case when (p = f and a = 3 is shown in Fig.^ and Figj^. From 
Figj^, in the absence of we note that the local extremes of N{p) approximately are regularly 
distributed at small times Gt ~ IOtt with the different values of 9. After that, the values of 
minimum for N{p) are increased when the time increases. Also, the time evolution curve of 
N(p) has a local extreme at the revival time tji. However, we note that, the value of N(p) is 
oscillated around its maxima for the different values of 9 and Gt. If we take into account the 
dispersive reservoir parameter, the amplitudes of the fast oscillations of the negativity measure 
diminish, and N{p) arrives to its minima (maxima) value at half-revival time. The negativity 
has a very strong sensitivity to ^ parameter. Therefore, when we take ^ = 0.1, see Figj^, the 
N(p) quite vanishes at particular times for the different values of 9. In this case, the dispersive 
reservoir at these particular times, the entanglement between the vibrational mode and qubit 
is completely destroyed. If N{p) = 0, the final state shown in Eqj^ disentangles completely 
and abruptly in just a finite time and evolves to the vacuum state, and its coherence is lost 
completely. 

To visualize the influence of detuning parameter and ^ parameter on the temporal evolution 
for negativity, the behavior of the negativity in off-resonance case with 0 = y, a = 3 and (p = f 
is displayed in Figj^ and Figgl. One can find that, the detuning parameter ^ 7 ^ 0 affects 
the amplitudes of fluctuations in the evolution of N{p) either by increasing or by decreasing 
on the case in which ^ = 0 for the different values of 9. Therefore, if 6 * = f, see Eigj^, the 
amplitudes of oscillations of tV(p) are less than their counterparts of the case ^ = 0. Because 
of the detuning parameter, some of oscillations are deteriorated, the local extremes of the 
negativity are pronounced and still apparent in this surface. In addition, with ^ ^ 0, the 
oscillations of the local extreme of the tV(p) are disappeared completely. When the dispersive 
reservoir is increased, say ^ = 0.1, see Eigj^, the amplitude of most oscillations are suppressed 
and tV(p) quite vanishes at particular times. This means that, after a particular time, the 
dispersive reservoir destroys the entanglement of the global total system state. Therefore, one 


11 







Figure 4: The time evolution of the entanglement N(p) in the resonance case for all 9 and 
off-resonance case for 9 = y, with ip = j and a = 3, for the different values of dispersive 
reservoir: ^ = 0 in (a, c) and ^ = 0.1 in (b, d). 
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can determine a particular region in which there is no entanglement between the qubit and the 
torsional resonator due to the complete destruction of coherence in the presence of 


4 Conclusion 

In this article, the system of SC phase qubit coupled to a torsional resonator is governed by a 
master equation. An analytical solution for the master equation is obtained with the general 
initial state. Several physical phenomena have been studied, such as qubit inversion, purity 
and negativity in the resonance and off-resonance cases for many different initial states. It is 
found that, the qubit inversion and purity, for the resonance and off-resonance cases, are quite 
sensitive to any change ^ parameter. By using the negativity the amount of the entanglement 
is studied. The dispersive reservoir destroys the entanglement of the global vibrational mode- 
phase qubit state. In addition, we found that in off-resonance case, the detuning parameter 
improve the influence of the damping on qubit inversion, purity and negativity. Therefore, one 
can determine a particular region in which there is no entanglement between the qubit and the 
torsional resonator due to the presence of T. 
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